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An integro-difTerential equation describing the angular distribution of beams is analyzed for a 
medium with random inhomogeneities. Beams are trapped because inhomogeneities give rise to wave 
localization at random locations and random times. The expressions obtained for the mean square 
deviation from the initial direction of beam propagation generalize the "3/2 law". 
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Fluctuations of the beam propagation direction in a randomly inhomogeneous medium 
are frequently observed in nature. Examples include random refraction of radio waves in 
the ionosphere and solar corona, stellar scintillation due to atmospheric inhomogeneities, 
and other phenomena. The propagation of a beam (of light, radio waves, or sound) in 
such media can be described as a normal diffusion process [H [2]. One extraordinary 
property predicted for random media - and later revealed - is the Anderson localization [3J . 
It brings normal diffusion to a complete halt. In the context of wave propagation in a 
random medium, the Anderson localization is caused by interference of waves resulting 
from multiple scattering [4j. When two waves propagating in opposite directions along a 
closed path are in phase, the resultant wave is more likely to return to the starting point 
than propagate in other directions. The properties of a randomly inhomogeneous medium 
vary not only from point to point, but also with time. Consequently, localization may 
take place both at random locations and at random times. Random localization affects 
diffusive light propagation in a random medium. An approach to description of this effect 
is developed in this paper. 

Suppose that the medium is statistically homogeneous and isotropic. Then, a beam 
propagating through the medium is deflected at random. Localization implies that the 
beam is trapped in some region. Since the trapped beam returns to the point where it 
was trapped, its propagation is "frozen" for some time. After that, a randomly deflected 
beam leaves the legion and propagates further until it is trapped in another region (or 
at a point) and the localization cycle repeats. The randomly winding beam path due to 
inhomogeneities is responsible for the random refraction analyzed in this study. 
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The angle 6 of deviation of a beam from its initial direction is characterized by a 
probability density Wa{0., cr), where a is the path travelled by the beam. Let us derive an 
integro-differential equation for the probability density. In contrast to rotational Brownian 
motion, the random walks analyzed here consist of random angle jumps A9i at points 
separated by segments of random length AcXj. Exact knowledge of the distributions of 
these random variables is not required. It is sufficient to assume that the angle jumps are 
independent random variables belonging the domain of attraction of Gaussian probability 
distribution. The random segment lengths Auj are also identically distributed independent 
random variables, with stable distribution is characterized by an index a. Since Acij is 
a nonnegative quantity, this distribution is totally asymmetric, and < a < 1. Recall 
that a random variable characterized by a probability distribution /(x), of this kind is 
described by the Laplace transform 



{s)= I exp{— sx} (i/(x) = exp{ — (Ax)"}, 
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where x > and A > [5]. The total path length is the sum of all Ao",. Both Acjj and 
A6'j are Markov processes. However, since the former is the master process with respect to 
the latter, the resultant process may not preserve the Markov property On account of 
convergence of distributions we can definitely pass from the discrete model to a continuous 
limit [7]. This leads to the diffusion equation 

w^{e^ a) - w^{e^ 0) = r ^.7-^^ f^i^^^^^^#^l - ^r-' da' , 
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where D is a diffusion coefficient and r(x) is the gamma function. The solution to this 
equation can be expressed as an integral transform of the probability density associated 
with rotational Brownian motion: 

/■oo 

Wo,i9, a) = / F^{z) Wi{9, a^z) dz , 
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where 



k=o ^'r(l — a — ka) 
The diffusion equation yields the mean 



cos^ = EJ-2Da") 



where 
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n=0 



is the Mittag-Leffier function. At large a, all beam directions are equiprobable. However, 
in contrast to normal diffusion (a = 1), a beam has to travel a longer path cr to reach 
this state. This process is somewhat analogous to "superslow" relaxation. 

Following the method developed in [8j, one can find the mean square of the distance r 
from the starting point to the observation point reached by the beam that has travelled 
an intricate path of length a through the medium: 

r^= ?^ ^(l- EJ-2Da'')). (1) 
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If Dcr -C 1, then 

^ ^ 2a- ^-^^] (2) 

^r(2a + l) r(3a + l)y ■ ^ ' 

If the 2 axis of a polar coordinate system is aligned with the initial beam direction, then 
the mean square of the distance passed by the beam along this axis is given by the formula 



1 

3D 



r(a + l) 



If Z^cr is small, then 

^ o 2a r 1 QDa' 



(3) 



(4) 



_r(2a + l) r(3a + l 

Now, the mean square deviation of the beam from its initial direction can be calculated 
by combining ([1]) with ([3]): 

'^ = 7^-^= - -J—(i-EJ-2Da'')) +^—(l-EJ-6Da'')). (5) 

If Da is small, then a generalized "3/2 law" is obtained [8]: 

2^2 



r(3a + r 



D^/V^^/l (6) 

The mean squares given by ([I]), ([3]) and increase as a" at large a. The case of a = 1 
corresponds to normal diffusion without wave localization. Thus, the approach developed 
here embraces classical results of the theory of beam propagation medium in a randomly 
inhomogeneous medium jTl |2] . 

Finally, it should be recalled that experimental deviations from the "3/2 law" (more 
precisely, from an exponent of 3/2 in the classical power law) were mentioned in |9]. 
However, they were attributed to systematic measurement errors, probably because of 
the lack of plausible interpretation. This problem can be revisited in view of the results 
obtained in this study. Moreover, new accurate experimental studies of beam propagation 
in appropriate randomly inhomogeneous media would be extremely useful for verifying 
the model proposed here. 
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